Neutrino Breakup of A = 3 Nuclei in Supernovae 
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We extend the virial equation of state to include 3 H and 3 He nuclei, and predict significant 
mass-three fractions near the neutrinosphere in supernovae. While alpha particles are often more 
abundant, we demonstrate that energy transfer cross-sections for muon and tau neutrinos at low 
densities are dominated by breakup of the loosely-bound 3 H and 3 He nuclei. The virial coefficients 
involving A = 3 nuclei are calculated directly from the corresponding nucleon- 3 H and nucleon- 3 He 
scattering phase shifts. For the neutral-current inelastic cross-sections and the energy transfer cross 
sections, we perform ab-initio calculations based on microscopic two- and three-nucleon interactions 
and meson-exchange currents. 
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I. INTRODUCTION 

Core-collapse supernovae (SN) are giant explosions of 
massive stars that radiate 99% of their energy in neutri- 
nos. Therefore, the dynamics and the neutrino signals 
can be sensitive to the details of neutrino interactions 
with nucleonic matter. At present, most SN simulations 
with detailed neutrino microphysics do not explode, but 
they may be close to successful explosions (for a status 
report, see [2, Q). Supernovae radiate electron, muon 
and tau neutrinos. Electron neutrinos can exchange en- 
ergy with matter via charged-current interactions. En- 
ergy transfer from muon or tau neutrinos, hereafter v x , 
is more difficult Q because neutrino-electron scattering 
has a small cross section, and neutrino-nucleon elastic 
scattering involves only a small energy transfer. 

Haxton and Bruenn proposed that v x can exchange 
energy via inelastic excitations of He and heavier nu- 
clei j4j, whereas Hannestad and Raffelt investigated the 
exchange of energy between v x and two interacting nuclc- 
ons i/ x NN — > i/ x NN [5[. Recently, Juodagalvis et al. have 
calculated detailed v x neutral-current cross-sections for 
A = 50 — 65 nuclei @, and Gazit and Barnea have pre- 
sented microscopic results for 4 He cross sections @, H|. 
These inelastic excitations can aid the transfer of neu- 
trino energy to the SN shock and can keep v x in thermal 
equilibrium to lower densities, resulting in the radiation 
of a lower energy v x spectrum First studies on the 
role of 4 He excitation for the shock revival were carried 
out by Ohnishi et al. (Ifjj . 

To evaluate the role of v x inelastic scattering one needs 
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both cross sections and detailed information on the com- 
position and other thermodynamic properties of nucle- 
onic matter. Models that describe the system with only 
a single average nuclear species may miss the contribution 
of less abundant nuclei with large cross sections. More- 
over, there are many fundamental connections between 
the equation of state and neutrino interactions. 

Nuclear statistical equilibrium (NSE) models predict 
abundances based on binding energies and the quantum 
numbers of nuclei. However, NSE models only treat 
approximately (or neglect) strong interactions between 
nuclei, and consequently break down as the density in- 
creases. We have recently developed a description of low- 
density nuclear matter (composed of neutrons, protons 
and alpha particles) in thermal equilibrium based on the 
virial expansion [III [l2| . The virial equation of state sys- 
tematically takes into account contributions from bound 
nuclei and the scattering continuum, and thus provides 
a framework to include strong-interaction corrections to 
NSE models. The virial equation of state makes model- 
independent predictions for the conditions [HI near the 
neutrinosphere, for low densities p ~ 10 11-12 g/cm 3 and 
high temperatures T ~ 4MeV. In particular, the result- 
ing alpha particle concentration differs from all equations 
of state currently used in SN simulations, and the pre- 
dicted large symmetry energy at low densities has been 
confirmed in near Fermi-energy heavy-ion collisions [l4l |. 
In addition, the long-wavelength neutrino response of 
low-density matter can be calculated consistently from 
the virial expansion (lEI ] . 

In this paper, we extend the virial expansion to in- 
clude 3 H and 3 He nuclei, and predict that the mass- 
three fraction can be significant (up to 10%) near the 
neutrinosphere. The second virial coefficients involving 
A = 3 nuclei are calculated directly from the correspond- 
ing nucleon- 3 H and nucleon- 3 He scattering phase shifts. 
While alpha particles are often more abundant due to the 
large binding energy (£4 = 28.3 MeV compared to E 3 ~ 
8 MeV) , we show that mass-three nuclei are important 
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for energy transfer, in particular for the more energetic 
muon and tau neutrinos with E Vx ~ 20 MeV. For neutri- 
nos with these energies, we find that the neutral-current 
inelastic energy transfer cross-sections and the neutrino 
energy loss for T > 4 MeV are dominated by the breakup 
of the loosely-bound 3 H and 3 He nuclei. Our predictions 
for the neutral-current inclusive inelastic cross-sections 
on mass-three nuclei are based on microscopic two- and 
three-nucleon interactions and meson-exchange currents, 
and include full final-state interactions via the Lorentz 
integral transform (LIT) method [l6| . 

This paper is organized as follows. In Section [TTJ we 
generalize the virial equation of state to include A = 3 
nuclei, and present results for the composition of low- 
density nuclear matter for various temperatures, densi- 
ties and proton fractions. In Section fill Al we calculate 
the inelastic 3 H and 3 He neutral-current cross-sections 
and energy transfer cross-sections. We combine our re- 
sults in Section IIIIBI and study the neutrino energy loss 
for conditions near the neutrinosphcrc. Finally, we con- 
clude in Section ITVl 



II. 



COMPOSITION OF LOW-DENSITY 
NUCLEAR MATTER 



In this section, we discuss the virial equation of state 
and present results for the composition of low-density 
nuclear matter including A = 3 nuclei (for details and 
additional thermodynamic results, see [17j|). 

To determine the abundance of A = 3 nuclei near the 
neutrinosphere in supernovae, we extend the virial ap- 
proach of Refs. [ll|, to include 3 H and 3 He nuclei. 
In the corresponding virial expansion, neutrons, protons, 
a particles, 3 H and 3 He nuclei are explicitly included, 
deuterons are included as a bound state contribution to 
the proton- neutron virial coefficient. We will explicitly 
consider deuterons and neutrino-deuteron scattering in 
future work. The equation of state is determined through 
an expansion of the pressure P in the fugacities (see for 
instance [18]) up to second order, 

P 2 
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where T is the temperature, Zi = e^ i+Ei ^ T is the fu- 
gacity (with chemical potential ^ and binding energy 
Ei), X{ = /rriiT is the thermal wavelength of parti- 
cle i, and b^ are the second virial coefficients describing 
interactions between particles i,j (bi = bo). We have 
calculated these virial coefficients from phases shifts and 
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TABLE I: Nucleon- 3 H and nucleon- 3 He virial coefficients for 
temperatures of T = 1 MeV to 10 MeV. 

low-energy scattering lengths, see Ref. [TtJ for details. 
The second virial coefficients b n , b pn , b a and b an are tab- 



ulated in [111 , and the additional virial coefficients involv- 
ing A = 3 nuclei are given in Table Q] from T = 1 MeV 
to 10 MeV. Due to a lack of p- 3 H scattering data, we as- 
sume that 6 p 3 H w 6„3 He . The effects of interactions with 
A = 3 nuclei are expected to be small at low densities (see 
also the hierarchy observed in Ref. [ll|) and therefore we 
neglect &3 He 3 He , &3 H 3 H , &3 Hc 3 H , 6 Q 3 Hc and 6 Q 3 H . 

The particle densities can be obtained by partial 
derivatives of the pressure with respect to the corre- 
sponding fugacities, 



d P 
dz~T 

By assuming chemical equilibrium, 

fi a 2/i n -(- 2/ip , 
fJ-m = 2^„ + Hp , 

/ i3 Hc = Hn + 2fl p , 



(2) 

(3) 
(4) 
(5) 



the corresponding A = 3,4 fugacities are determined by 
the proton and neutron fugacities and the correspond- 
ing binding energy: Z3 H = z^z p e E3 v-f T for example. For 
a given baryon density nt, and proton fraction Y p , the 
fugacities z p and z n are then determined implicitly from 

rib = n p + n n + An a + 3m Hc + 3ns H , (6) 
Y p = (rip + 2n a + 2ri3 Hc + n3 H )/n b . (7) 

The resulting mass fractions (x^ — AiUi/rib) are shown 
in Fig. [1] for neutrinosphere densities and temperatures, 
and various proton fractions. 

For a density of 10 12 g/cm 3 , we find that the total 
mass-three fraction is significant, up to a;3 H + x3 Hc 1=3 0.1 
in symmetric matter. Moreover, for temperatures T > 
5 MeV (and increasingly important at low proton frac- 
tions), the 3 H mass fraction is larger than the a particle 
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FIG. 1: (Color online) Mass fractions of nucleons and A — 3, 4 nuclei in chemical equilibrium as a function of temperature 
T. The top and bottom rows correspond to a density of 10 11 g/cm 3 and 10 12 g/cm 3 respectively, and from left to right the 
proton fractions are Y p = 0.1, 0.3 and 0.5. For this temperature range, at 10 11 g/cm 3 the neutron fugacity is z n < 0.1, 0.05, 0.01 
for Y p = 0.1,0.3 and 0.5 respectively, all other fugacities are < 0.01 for all shown proton fractions. At the higher density, 
10 12 g/cm 3 , the neutron fugacity is z n < 0.7, 0.35, 0.04 for Y p = 0.1, 0.3 and 0.5 respectively, all other fugacities are < 0.06 for 
all shown proton fractions. 



fraction. We also observe that at lower temperatures, 
where many of the protons are bound in alpha particles 
(and heavy nuclei), 3 H nuclei are more abundant than 
free protons. Here, the rate of electron capture may be 
dominated by capture on 3 H, since the capture on 4 He 
has low cross sections and there are few free protons. 
The contribution of mass-three nuclei to charged-current 
interactions will be left to future work [11]. Finally, we 
emphasize that heavier nuclei and larger clusters become 
important as the a particle fraction saturates at very low 
temperature. 



III. NEUTRINO- 3 H AND - 3 He BREAKUP: 
NEUTRAL-CURRENT CROSS-SECTIONS 

The calculation of the neutral-current inclusive in- 
elastic cross-sections on 3 H and 3 He follows Ref. Q. 
We solve the three-nucieon problem based on the Ar- 
gonne Vw nucleon-nucleon [20| and the Urbana IX three- 
nucleon [21j interactions. Neutrino scattering on A = 3 
nuclei only induces transitions to continuum states, since 
3 H and 3 He have no excited states. Hence, a correct 
description must include breakup channels and final- 
state interactions among the three nucleons. We in- 
clude these via the LIT method [l6[, which uses an 
integral transform with a lorentzian kernel to reduce 



the continuum problem to a bound-statc-like problem. 
The resulting Schrodinger-like equations are solved us- 
ing the effective interaction hyperspherical harmonics 
(EIHH) approach [22|, [H| . The combination of these ap- 
proaches converges rapidly and yields a numerical preci- 
sion of less than a percent for few-body reaction cross- 
sections @, H, [H], [23 . The energy transfer due to elastic 
scattering is low, u> ~ T 2 /m, and therefore we include 
only breakup channels in our calculations. 

Since the energy scale of SN neutrinos is much smaller 
than the mass of the iv-boson, the neutrino-nucleus in- 
teraction can be approximated by an effective current- 
current Hamiltonian. The neutrino current is straight- 
forward and results in kinematical factors to the cross 
section. The standard model dictates only the formal 
structure of the nuclear neutral-current: 

J° = (1 - 2 sin 2 0w) y < + y < - sin 2 6 W J% , (8) 

where the superscripts A, V denote axial and vector cur- 
rents. For supernova neutrinos, chiral effective field the- 
ory (EFT) of nucleons and pions offers a consistent ap- 
proach to nuclear interactions and electro- weak currents. 
For historical reasons, the present calculation uses con- 
ventional two- and three-nucleon interactions and EFT 
currents, but future applications can be fully based on 
chiral EFT. The current approach has been applied to 
study electro- weak reactions on A = 2, 3, 4 nuclei [1, Hf|- 
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We use chiral EFT meson-exchange currents (MEC) 
at next-to-next-to-next-to-leading order. The MEC are 
based on a momentum expansion in Q/A, where Q ~ 
10 — 20 MeV is the typical energy in our process of inter- 
est, and the cutoff A is of the order of the EFT breakdown 
scale. Here, we follow Park et al. [IB] and vary the cut- 
off over the range A = 400 — 800 MeV. In configuration 
space, the MEC are obtained from a Fourier transform 
of propagators with a cutoff A. This leads to a cutoff de- 
pendence, which is renormalized by a cutoff-dependent 
countcrtcrm. In the present case, all other low-energy 
coefficients can be determined from pion-nucleon scatter- 
ing. The counterterm d r (A) characterizes the strength of 
a two-nucleon contact operator and has been matched to 
the triton half-life over this cutoff range. As a check, we 
reproduce the cutoff dependence d r (A) of Ref. (2r| . 



A. Inelastic cross-sections and energy transfer 

The calculated cross sections are averaged over energy 
and angle, assuming a Fcrmi-Dirac distribution for the 
neutrinos with zero chemical potential, temperature T„, 
and neutrino momentum k. 



f{T V) k) = 



N 



T„ 3 e k / T * + 1 



(9) 



where ./V -1 = 2 YlnLi{~ l) n+1 /^ 3 is a normaliza- 
tion factor. The quantities of interest are the 
temperature-averaged cross-sections and energy transfer 
cross-sections: 



( a )T„ — [ dej f dh f(T v ,ki) 

Jt*.* J dk f 



dk f 

f°° f do 
{ua) Tu = I du / dki f{T v ,ki)uj—— , (11) 

•W J dk f 



(10) 



where kij are the initial and final neutrino energy, 
uj = ki — kf is the energy transfer, and u>th denotes the 
threshold energy of the breakup reaction. In Table m 
we present results for the averaged neutrino and anti- 
neutrino neutral-current inclusive inelastic cross-sections 
and energy transfer cross-sections as a function of the 
neutrino temperature. The difference between the 3 H 
and the 3 He cross sections reflects the difference in 
thresholds between the two nuclei. The mirror symme- 
try between both nuclei is restored with higher neutrino 
energy. The leading contributions to the cross section 
are the axial E^, Mj 4 and E£ multipoles. The relative 
importance of these multipoles varies as a function of the 
momentum transfer, and thus as a function of the neu- 
trino temperature. In comparison to inelastic excitations 
of 4 He studied in Ref. § , we find that the cross sections 
are about a factor 20 and 10 times larger at temperatures 
4 MeV and 6 MeV respectively (for the mean values of 3 H 
and 3 He), and the energy transfer cross-sections are 8 and 
2 times larger. 
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TABLE H: Averaged neutrino- and anti-neutrino- 3 H and - 3 He 
neutral-current inclusive inelastic cross-sections per nucleon 
(A=3), {o)t u = (&v + °""^)t„ (left columns), and energy 
transfer cross-sections, (uxj)t v = jji^* 7 " + W0y}T„ (right 
columns), as a function of neutrino temperature T v , in units 
of 10~ 42 cm 2 and 10~ 42 MeVcm 2 respectively. 



At low-momentum transfer, the Gamow- Teller oper- 
ator dominates for the cross section, consequently the 
MEC have a large effect of about 16% at a tempera- 
ture of 1 MeV. At higher-momentum transfer, higher- 
order multipoles start to play an important role. Due to 
spatial symmetry, the MEC contribution to these mul- 
tipoles is small and the overall effect of MEC decreases 
rapidly to < 2% for temperatures above 4 MeV. While 
not directly important here, the asymmetry between the 
scattering of neutrinos and anti-neutrinos increases with 
temperature: The difference in the energy transfer grows 
gradually from 3% for a neutrino temperature of 3 MeV 
to > 50% for 10 MeV temperatures. Finally, the cutoff 
dependence of these observables is < 2% for 1 MeV and 
< 1% for higher temperatures. This validates our calcu- 
lations. We thus estimate the precision of the predicted 
cross sections to be a few percent, which also includes 
estimates of the numerical accuracy. 



B. Neutrino energy loss due to inelastic scattering 

We can combine the energy transfer cross-sections with 
the A — 3,4 mass fractions of Section |TT] to calculate 
the neutrino energy loss due to inelastic excitations of 
A = 3,4 nuclei. The neutral-current cross-sections on 
4 He are taken from Ref. [|[, which is based on the same 
microscopic input. A neutrino of energy E v will lose 
energy to inelastic excitations, and heat the matter, at a 
rate dE^/dx given by 



dE v 
dx 



= 3 H. 3 Hc, 4 Hc 



(12) 



To explore the effect of mass-three nuclei on the energy 
loss, we assume the neutrino energies are characterized 
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FIG. 2: (Color online) Neutrino energy loss dE v /dx for inelastic excitations of A = 3, 4 nuclei as a function of the matter 
temperature T at a density of 10 12 g/cm 3 . We assume that the the neutrino energies are characterized by a Fermi-Dirac 
distribution with a temperature T„ = 6 MeV. The contributions from 3 H, 3 He and 4 He nuclei, and the total neutrino energy 
loss are shown for proton fractions Y p = 0.1, 0.3 and 0.5. 



by a Fermi-Dirac spectrum of temperature T„, while the 
low-density matter may have a lower temperature T. For 
simplicity, we neglect the energy transfer from nuclei to 
neutrinos required by detailed balance. This is strictly 
correct only in the limit T <C T„. 

In Fig. [21 the neutrino energy loss due to inelastic 
scattering is shown for a density of 10 12 g/cm 3 and neu- 
trino temperature T„ = 6 MeV, as a function of the 
matter temperature for various proton fractions. For 
T > 4 MeV, the energy loss is dominated by the con- 
tributions from 3 H nuclei. The total abundance of A = 3 
nuclei depends only weakly on the proton fraction (see 
Fig. [1]) , which is reflected in the weak dependence of the 
neutrino energy loss as a function of proton fraction. Fi- 
nally, for lower densities, mass-three nuclei are less abun- 
dant (see Fig. [TJ, and therefore also their contributions 
to the neutrino energy loss. 



IV. CONCLUSIONS 



free protons at low temperatures, since many of the pro- 
tons are bound in alpha particles (and heavy nuclei). In 
these regions, the rate of electron capture may be domi- 
nated by capture on 3 H nuclei [l9l ]. 

While alpha particles are often more abundant, we 
have shown that the loosely-bound 3 H and 3 He nuclei 
dominate the energy transfer at low densities through in- 
elastic excitations, and are therefore especially important 
for energy transfer from muon and tau neutrinos. Our 
new results for neutrino- 3 H and - 3 He neutral-current in- 
clusive inelastic cross-sections and energy transfer cross- 
sections are based on microscopic two- and three-nucleon 
interactions and meson-exchange currents, which repro- 
duce the triton half-life. All breakup channels and 
full final-state interactions were included via the LIT 
method. For temperatures T ~ 4 MeV, the predicted 
energy transfer cross-sections on mass-three nuclei are 
approximately one order of magnitude larger compared 
to inelastic excitations of 4 He nuclei. 



The virial expansion provides a systematic approach 
to low-density nuclear matter in thermal equilibrium, 
in particular for the conditions near the neutrinosphere 
with densities p ~ 10 11-12 g/cm 3 and temperatures T ~ 
4 MeV. In this paper, we have extended the virial equa- 
tion of state of Refs. [ll],[l2| (for matter composed of neu- 
trons, protons and alpha particles) to include 3 H and 3 He 
nuclei. We have made model-independent predictions for 
the abundance of 3 H and 3 He nuclei and predicted that 
their mass fractions can be significant near the neutri- 
nosphere. Our results are directly based on nucleon- 3 H 
and nucleon- 3 He scattering phase shifts. In addition, it 
is interesting that 3 H nuclei can be more abundant than 



Using the virial abundances and the microscopic en- 
ergy tranfer cross-sections, we have found that mass- 
three nuclei contribute significantly to the neutrino en- 
ergy loss due to inelastic excitations for T > 4 MeV. In- 
elastic excitations of 3 H and 3 He nuclei can be more im- 
portant than He nuclei at high temperatures, low pro- 
ton fraction or higher densities. To fully assess the role 
of neutrino breakup of A = 3 nuclei, our predicted abun- 
dances and neutral-current cross-sections should be in- 
cluded in SN simulations. The model independence of 
the virial equation of state and the accuracy of the pre- 
dicted cross sections can help to improve the theoretical 
nuclear microphysics for SN simulations. 
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